Observation of a Pree-Shercliff-Layer Instability in Cylindrical Geometry 
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We report on observations of a free-Shercliff-layer instability in a Taylor-Couette experiment using 
a liquid metal over a wide range of Reynolds numbers, Re ~ 10 3 — 10 6 . The free Shercliff layer is 
formed by imposing a sufficiently strong axial magnetic field across a pair of differentially rotating 
axial endcap rings. This layer is destabilized by a hydrodynamic Kelvin-Helmholtz-type instability, 
characterized by velocity fluctuations in the r — 9 plane. The instability appears with an Elsasser 
number above unity, and saturates with an azimuthal mode number m which increases with the 
Elsasser number. Measurements of the structure agree well with 2D global linear mode analyses 
and 3D global nonlinear simulations. These observations have implications for a range of rotating 
MHD systems in which similar shear layers may be produced. 



The destabilization of rotating sheared flows by an ap- 
plied magnetic field in magnetohydrodynamics (MHD) 
is a topic with astrophysical and geophysical implica- 
tions, and has been the subject of a number of exper- 
imental and theoretical efforts. Such destabilization can 
be caused by the magnetorotational instability (MRI), in 
which a magnetic field of sufficient amplitude can desta- 
bilize Rayleigh-stable rotating sheared flows In this 
Letter, we report the observation of an instability which, 
like the MRI, appears in a sheared rotating fluid when a 
magnetic field is applied. But rather than playing a role 
in the dynamics of the instability, as in the case of the 
MRI, the magnetic field here acts to establish free shear 
layers which extend from axial boundaries and which are 
subject to a hydrodynamic instability. 

Hartmann and Shercliff laid the groundwork in under- 
standing the effect of magnetic fields on shear layers in 
conducting fluids. Hartmann studied boundary layers 
normal to an external applied field , and Shercliff ex- 
tended his analysis to include boundary layers parallel to 
the applied field Free Shercliff layers can be estab- 
lished in rotating MHD systems when the line-tying force 
of an axial magnetic field extends a discontinuity in an- 
gular velocity at an axial boundary into the bulk of the 
fluid. These shear layers are similar to the Stewartson 
layers that extend from discontinuous shearing bound- 
aries in rapidly rotating systems [4( , but for the free Sher- 
cliff layer discussed here, it is the magnetic field tension 
rather than the Coriolis force that leads to equalization 
of the angular velocity in the axial direction. 

Free Shercliff layers were first realized experimentally 
by Lehnert in a cylindrical apparatus with a free sur- 
face at the top and a rotating ring at the bottom axial 
boundary [5( . Lehnert observed the formation of vortices 
at the location of the shear layers, though he attributed 
their formation to discontinuities in the free surface at the 
shear layer location rather than to the shear itself. These 
layers were then described analytically by Stewartson [f| 
and Braginskii [3] • The formation of free Shercliff layers 



in a cylindrical Taylor-Couette device has been predicted 
computationally 8], but these simulations were axisym- 
metric and thus incapable of evaluating the stability of 
these shear layers to nonaxisymmetric perturbations. 

Both free Shercliff layers and Stewartson layers can be 
present at the tangent cylinder of spherical Couette sys- 
tems. The Kelvin-Hclmholtz destabilization of these lay- 
ers has been studied extensively through computation [9|- 
Ulj |. Stewartson layers have been observed experimen- 
tally in spherical and cylindrical geometry and are found 
to be unstable to nonaxisymmetric modes, which is con- 
sistent with simulations EH, EH . 

The Princeton MRI experiment is a Taylor-Couette ap- 
paratus consisting of two coaxial stainless steel cylinders 
as shown in Fig.Q] The gap between the cylinders is filled 
with a GalnSn eutectic alloy which is liquid at room tem- 
perature. Differential rotation of the cylinders sets up a 
sheared rotation profile in the fluid. If the cylinders were 
infinitely long, the fluid between the cylinders would as- 
sume an angular velocity VL at a radius r matching the 
ideal Couette solution in steady state, f2(r) = a + b/r 2 . 
The constants a and b are found by matching the solu- 
tion to the imposed rotation rates at the inner and outer 
cylinder boundaries. In conventional Taylor-Couette de- 
vices, the endcaps are typically corotatcd either with the 
inner or outer cylinder. This produces strong secondary 
circulation and angular momentum transport to the ax- 
ial boundaries, resulting in a deviation from the ideal 
rotation profile 15]. A novelty of this apparatus is the 



configuration of the axial endcaps, each of which is split 
into two differentially-rotatable acrylic rings, giving four 
independent rotation rates: those of the inner cylinder, 
outer cylinder, inner rings, and outer rings. In previ- 
ous experiments using water as the working fluid, this 
configuration was very effective at reducing the influence 
of the axial boundaries, allowing the generation of quies- 
cent flows in the bulk of the fluid with Reynolds numbers 
Re = f2iri(r2 — r i)/v above 10 6 [16|. The experimental 
parameters are shown in Table H] 
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FIG. 1: Diagram of Princeton MRI experiment. Each endcap 
is split into an inner ring (IR) and an outer ring (OR). Differ- 
ential rotation of these rings produces a discontinuity in the 
angular velocity at the axial boundary. Overlaid on the right 
half of the figure is a plot of the shear (r/Q)(dQ/dr) from a 
nonlinear MHD simulation with differential rotation between 
the endcap rings and a strong axial magnetic field [14J. The 
free Shercliff layers are the regions of strong negative shear 
extending from the interface between the rings. 

Fluid velocities are measured with an ultrasound 
Doppler velocimetry (UDV) system[T3, EH- Ultrasonic 
transducers are mounted on the outer cylinder at the 
midplane of the experiment. A transducer aimed radially 
and others aimed tangential to the inner cylinder allow 
determination of the radial and azimuthal velocity com- 
ponents. Two tangential transducers aimed identically 
but separated azimuthally by 90° provide information 
about azimuthal mode structure. 

A set of six solenoidal coils applies an axial magnetic 
field to the rotating fluid. Fields below 800 Gauss can be 
applied indefinitely, while the application time for higher 
fields is limited by the resistive heating of the coils. An 
array of 72 magnetic pickup coils placed beyond the outer 
cylinder measures dB r /dt. 

Experiments were run using both Rayleigh- stable and 
-unstable flow states. The Rayleigh- stable states had 
component rotation speeds in the ratio [1.0, 0.55, 0.1325, 
0.1325] for the inner cylinder, inner ring, outer ring, 
and outer cylinder, respectively. The ideal Couette so- 
lution for these inner and outer cylinder speeds satis- 
fies Rayleigh's stability criterion that the specific angular 
momentum increase with radius: d(r 2 Q)/dr > 0. The 



TABLE I: Parameters of the apparatus [19J and liquid metal 
working fluid [2p| ]. 



Parameter 


symbol 


value 


units 


Height 


h 


27.9 


cm 


Inner cylinder radius 


n 


7.06 


cm 


Outer cylinder radius 


rz 


20.3 


cm 


Density 


P 


6.36 


g/cm 3 


Kinematic viscosity 


V 


2.98 x 10~ 3 


cm 2 /s 


Magnetic diffusivity 


V 


2.57 x 10 3 


cm 2 /s 


Inner cylinder rotation rate 


Hi 


0.25 - 800 


rpm 


Axial magnetic field 


B 


- 4500 


Gauss 



ring speeds were chosen empirically to generate an az- 
imuthal rotation profile in the hydrodynamic case that 
closely matches the ideal Couette profile at the midplane. 
The Rayleigh- unstable states were generated using com- 
ponent speeds in the ratio [1.0, 1.0, 0, 0]. These flows vi- 
olate Rayleigh's criterion and exhibit large velocity fluc- 
tuations in the absence of a magnetic field. 

A single run of this experiment starts with an acceler- 
ation phase of two minutes, during which the sheared 
azimuthal flow develops. The axial magnetic field is 
then applied, initially resulting in the damping of hy- 
drodynamic fluctuations. If the magnetic field is strong 
enough to satisfy the requirement that the Elsasser num- 
ber A = B 2 /4Trpr]An, > 1, where Afi is the difference 
between the inner- and outer-ring rotation rates, the in- 
stability grows up as a large-scale coherent mode. It man- 
ifests itself as a fluctuation in the radial velocity and az- 
imuthal velocity, where significant perturbations of more 
than 10% of the inner cylinder speed are observed. An ul- 
trasonic transducer inserted on a probe and aimed axially 
at an endcap did not measure axial velocity fluctuations 
when the instability was excited, suggesting that the flow 
due to the instability is mainly in the r — 9 plane. Cor- 
related magnetic fluctuations are observed at the highest 
rotation rates and applied fields. The instability develops 
on both the Rayleigh-stable and -unstable backgrounds, 
and typical mode rotation rates exceeds the outer cylin- 
der rotation rate 0,2 by ~ 0.1(Sli — f^)- 

The instability was observed over a range of more than 
3 orders of magnitude in rotation rate in the Rayleigh- 
unstablc configuration, as shown in Fig. [2] with Re — 
820—2.6 x 10 6 . The instability is present even with a mag- 
netic Reynolds number Rra = Oifi(r2 — fij/v ~ 10~ 3 , 
indicating an inductionless mechanism in which induced 
magnetic fields are dynamically unimportant. 

For A of order one, the primary azimuthal mode num- 
ber at saturation is m = 1, with phase- locked higher- 
order mode numbers typically present at a smaller am- 
plitude. The measured mode structure is shown in Fig. [3] 
It is common for an to = 2 mode to grow up before an 
to = 1 dominates at saturation. High-A scenarios at very 
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FIG. 2: Stability diagram for the Rayleigh-unstable back- 
ground flow state. The area of the circles is proportional to 
the power in the dominant Fourier harmonic measured by a 
tangential transducer at r = 19.2cm, normalized to the square 
of the inner cylinder speed. The 'x's indicate stability. The 
stability diagram for the case starting from a Rayleigh-stable 
background is similar, but was measured over a smaller range 
of speeds. The inset plot shows a sample time trace of the 
velocity measured at one point in the flow, with the magnetic 
field applied in the region between the dashed vertical lines. 
(tti =200 rpm, B =2900 G.) 

slow rotation rates show that m at saturation increases 
as A increases. Primary mode numbers up to m = 5 
have been observed with A = 127 at a rotation rate of 
0.25 rpm. 

The necessity of shear at the axial boundary has been 
verified experimentally. Experiments were performed 
with the components rotating in the standard Rayleigh- 
stable configuration, but with a number of different in- 
ner ring speeds. The critical magnetic field for instability 
varied with the differential rotation between the endcap 
rings as expected. When the inner rings and outer rings 
corotated, the instability was not observed. 

The free shear layer has been measured experimentally 
at low Re and high A where it penetrates to the midplane 
of the experiment as shown in Fig. |4j The width of the 
layer measured at a time just before the onset of insta- 
bility is consistent with the expected width scaling for 
a Shercliff layer 5 ~ 1/ vM, where the Hartmann num- 
ber M = Bl/y/^Trprjv and I = T2 — T\ is a characteristic 
length. The onset of the instability is associated with a 
decrease in the mean shear in this layer. 

Nonlinear numerical MHD simulations have been per- 
formed with the HERACLES code [211 modified to in- 
clude finite viscosity and resistivity [14j. The simula- 
tions were performed in the experimental geometry with 
a 200x64x400 grid in f, 6, and z, with Re = 4000 and a 
range of Rm and M. These simulations show the forma- 
tion of the free Shercliff layer extending from the disconti- 
nuity at the axial boundaries, as shown in Fig.[TJ The ax- 




FIG. 3: Comparison of measured unstable mode with results 
from simulation. All are contour plots of azimuthal veloc- 
ity at the midplane with the m = contribution subtracted. 
Red indicates positive velocity (counterclockwise), and blue 
indicates negative velocity (clockwise). Upper left: Experi- 
mental measurement with A = 1.4, reconstructed from pro- 
jecting the time behavior of the azimuthal velocity as mea- 
sured by one UDV transducer onto the r — 6 plane. Upper 
right: Experimental measurement with A = 50. Lower left: 
Growing m — 1 mode produced by a hydrodynamic linear sta- 
bility analysis of an axially independent shear layer. Lower 
right: Unstable mode from nonlinear MHD calculation with 
Re = 4000, Rm = 10, and A = 1. 

ial length of the shear layer scales with \ZA, which seems 
to arise from a competition of magnetic forces, which act 
to extend the shear layer into the fluid, and poloidal cir- 
culation generated by the axial boundaries, which acts 
to disrupt the free shear layer. The simulations also pro- 
duce an instability requiring A > 1 for onset and suggest 
that a minimum penetration depth of the shear layer is 
required for development of the instability. Like the ex- 
perimental observations, the unstable modes exhibit a 
spiral structure, and a cascade is observed from higher 
azimuthal mode number during the growth phase of the 
instability to a dominant m = 1 at saturation. 

A global linear stability analysis was performed to in- 
vestigate unstable modes in the experimental geometry. 
The analysis found eigenvalues of the linearized nonideal 
MHD equations discrctized across 2048 grid cells in the 
radial direction, assuming sinusoidal azimuthal depen- 
dence with a specified mode number and no axial depen- 
dence. Unstable hydrodynamic solutions were sought for 
realistic fluid parameters and for a zeroth order, back- 
ground rotation profile consisting of a free shear layer 
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FIG. 4: Angular velocity versus radius at the midplane for 
several values of applied fields at a rotation rate of 0.5 rpm 
in the Rayleigh-unstable configuration. The dashed line in- 
dicates the radial position of the split between the axial end- 
caps. The inset plot shows measurements of the shear layer 
thickness from experiments at 0.25 rpm (green) and 0.5 rpm 
(blue) versus the Hartmann number M. M = 55%/A at 0.5 
rpm, and M = 39VA at 0.25 rpm. The solid line indicates 
the expected Shercliff layer width scaling S ~ d/\/ r M, where 
the constant d = 62cm has been chosen to match the data. 



represented by a hyperbolic tangent centered between 
the inner and outer cylinders. Angular velocity profiles 
with a sufficiently narrow shear layer were found to be 
hydrodynamically unstable to nonaxisymmetric Kelvin- 
Helmholtz modes with a similar structure to those ob- 
served experimentally. The most unstable mode number 
increases with decreasing shear layer width, similar to 
the experimental observations of the saturated states. 

The results presented here describe a minimum mag- 
netic field required for onset of the instability. Simu- 
lations have shown that a sufficiently strong magnetic 
field will restabilize this instability, similar to the simu- 
lation results in spherical geometry Experimentally, 
the decreasing saturated amplitude with increasing field 
at small rotation rates, shown in Fig. [3J suggests that 
this critical field strength is being approached. But the 
limits on controllable slow rotation and on the availabil- 
ity of strong magnetic fields precluded verification of the 
complete restabilization. 

This free-Shercliff-layer instability exhibits strong sim- 
ilarities to the expected behavior of the standard MRI 
in a Taylor- Couette device because in both cases a mag- 
netic field acts to destabilize otherwise stable flow and in 
both cases the associated angular momentum transport 
results in a large modification to the azimuthal velocity 
profile. But this instability is a hydrodynamic instability 
on a background state established by the magnetic field 
and is present with Rm -C 1. While there are induction- 
less relatives of the standard MRI, such as the so-called 
HMRI which relies on azimuthal and axial applied mag- 



netic fields 22], |23j , the unimportance of induction here is 
in stark contrast to the requirement of a finite minimum 
Rm for the standard MRI in an axial magnetic field. 

These results have particular relevance to other MHD 
experiments in which similar shear layers may be estab- 
lished. A spherical Couette MHD experiment produced a 
nonaxisymmetric instability with a ppli ed magnetic field 
that was claimed to be the MRI [24J. However, sub- 
sequent simulations have attributed those observations 



to hydrodynamic instability of free shear layers [25 . 
similar to the observations that we report. We expect 
that other cylindrical devices, such as the PROMISE 2 
experiment [27J, could produce this instability. But the 
critical value of A will likely change for experiments with 
different geometric aspect ratios. 

The free-Shercliff-layer instability is not expected to 
impact the study of the MRI in this device since the 
magnetic fields required for the MRI are weaker than 
those required for the Shercliff layer instability at MRI- 
relevant speeds [14] . 
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